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Abstract 


Two  robust  adaptive  nonlinear  controller  designs  are  presented  in  this  report  for  control  of  an 
air-breathing  hypersonic  vehicle  in  the  cruise  phase  of  flight.  The  first  type  of  controller  uses 
dynamic  inversion  and  the  second  one  is  obtained  using  newly  developed  state-constrained 
generalized  dynamic  inversion  technique  with  the  help  of  Lyapunov  theory.  Furthermore, 
the  robustness  of  both  controllers  is  enhanced  by  augmenting  them  with  a  fast  disturbance 
observer. 

The  controller  is  derived  using  dynamic  inversion  technique,  by  transforming  nonli¬ 
near  system  dynamics  into  linear  system  dynamics  with  the  help  of  transformation  matrix. 
Further,  the  control  expression  is  obtained  using  stable  linear  error  dynamics,  which  ensures 
the  asymptotic  stability  of  the  error.  However,  a  perfectly  known  system  is  assumed  while 
designing  the  controller  using  dynamic  inversion,  which  is  very  difficult  to  achieve  in  the  case 
of  air-breathing  hypersonic  vehicle.  Hence,  constrained  neuro-adaptive  dynamic  inversion 
technique  with  Jacobian  learning  approach  is  proposed  in  this  report.  In  this  technique,  an 
unknown  disturbance  function  is  learned  by  the  controller  along  with  the  Jacobian  matrix  of 
the  unknown  disturbance  to  ensure  fewer  transients.  Further,  barrier  Lyapunov  function  is 
used  in  the  formulation  of  neuro-adaptive  dynamic  inversion  technique  to  constrain  the  error 
in  disturbance  learning. 

The  state-constrained  generalized  dynamic  inversion  technique  is  formulated  using 
Barrier  Lyapunov  function  in  this  report.  Further,  the  derived  controller  is  robust  enough 
to  maintain  the  states  within  the  constrained  bounds  in  presence  of  unknown  disturbances 
or  model  uncertainties.  However,  the  perfect  tracking  is  achieved  by  augmenting  the  state- 
constrained  generalized  dynamic  inversion  with  constrained  neuro-adaptive  Jacobian  matrix 
learning  based  disturbance  learning  mechanism. 
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Further,  six  degree-of-freedom  winged-cone  air-breathing  hypersonic  vehicle  model 
is  used  in  this  report  to  demonstrate  the  performance  of  derived  controllers.  Moreover,  the 
uncertainties  are  assumed  in  aerodynamic  coefficients  and  inertia  components  of  the  vehicle 
model.  Further,  it  has  been  shown  that  the  proposed  controllers  are  robust  enough  to  take 
care  of  model  uncertainties  or  any  bounded  unknown  external  disturbance  function  affecting 
the  air-breathing  hypersonic  vehicle  dynamics.  Along  with  the  stabilization  of  uncertain 
air-breathing  hypersonic  vehicle,  the  controllers  also  prevent  inlet  unstart.  It  has  also  been 
shown  that  the  proposed  controllers  can  perfectly  tolerate  up  to  30%  parameter  uncertainties 
and  disturbances.  Moreover,  using  Monte- Carlo  simulations,  the  quantitative  results  of  the 
performance  of  proposed  controllers  in  presence  of  model  uncertainties  are  also  presented 
in  this  report.  Finally,  the  applicability  of  the  proposed  control  techniques  is  illustrated  by 
promising  simulation  results. 
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Introduction 


1.1  Motivation 

Air-breathing  hypersonic  vehicles  (AHV)  offer  distinct  performance  advantages  over  rocket- 
based  systems  for  space  access  vehicles.  The  major  advantage  of  using  air-breathing  engine  is 
that  the  extra  oxidizer  is  not  required  for  propulsion  (rather  oxygen  from  the  atmosphere  is 
used  as  the  oxidizer).  Moreover,  there  are  no  moving  parts  in  the  engine  such  as  compressors 
and  turbines,  Hence,  the  propulsive  system  is  easily  integrated  with  the  vehicle.  Further,  the 
speed  of  the  hypersonic  vehicle  is  very  high  (more  than  five  times  of  sound  speed),  therefore, 
it  takes  very  less  time  to  reach  the  target/ destination.  However,  these  performance  advantages 
are  dependent  upon  advances  in  current  state-of-the-art  technologies  in  many  critical  areas, 
one  of  which  happens  to  be  the  control  system  design. 

The  control  and  stabilization  of  AHV  are  very  challenging  because  its  aerodynamic  cha¬ 
racteristics  are  different  from  the  conventional  (subsonic  or  supersonic)  air  vehicles.  Moreover, 
High  velocity  of  AHV  results  into  highly  sensitive  flight  conditions  which  result  in  instability. 
This  instability  occurs  mainly  due  to  the  constant  center  of  pressure  position  irrespective  of 
the  angle  of  attack  and  altitude.  The  high  velocity  of  the  vehicle  also  results  in  weakly  damped 
transient  oscillations  of  the  airframe  [Coleman  and  Faruqi,  2009]. 

Furthermore,  the  engine  performance  of  hypersonic  vehicles  is  quite  sensitive  to  the 
angle  of  attack  and  the  sideslip  angle.  Deviation  of  these  angles  can  even  lead  to  the  disastrous 
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condition  of  engine  unstart.  Moreover,  at  hypersonic  speed,  the  behavior  of  the  vehicle  is 
highly  nonlinear  and  coupled  [Schmidt,  1993].  The  high  temperature  and  boundary  layer 
effect  on  control  surfaces  place  constraint  on  control  surface  deflection  limits.  This  restricted 
control  deflection  limit  places  a  great  difficulty  to  a  controller  designer.  This  may  lead  to 
quick  control  saturation  as  well  as  less  aerodynamic  damping,  leading  to  additional  difficulties 
for  control  design  [McRuer,  1991,  Fidan  et  al.,  2003].  Furthermore,  hypersonic  speeds  cause 
the  so-called  “path-attitude  decoupling"  phenomenon,  which  is  the  resistance  of  the  high 
momentum  of  the  hypersonic  vehicle  to  the  change  in  the  desired  flight  path  [Mirmirani 
et  al.,  2005b].  Hence,  the  actual  flight  path  gets  fairly  decoupled  from  the  pitch  attitude  at 
hypersonic  speeds. 

Furthermore,  Static  stability  margins  decrease  as  the  Mach  number  increases  in  AHV. 
Hence,  in  order  to  have  reasonable  static  margins  at  hypersonic  speeds,  one  needs  to  accept 
high-static  stability  margins  [Mirmirani  et  al.,  2005b].  Moreover,  the  stability  of  the  AHV 
can  not  be  defined  only  using  static  margin  because  of  changes  in  center  of  gravity  position 
(due  to  fuel  burn-off).  Hence,  the  active  control  augmentation  is  required  to  ensure  the 
stability  of  the  vehicle.  The  controller  design  problem  is  further  complicated  because  of  the 
fact  that  hypersonic  aerodynamic  parameters,  as  predicted  from  ground  tests  or  theoretical 
computational  methods,  do  not  reflect  the  actual  flight  parameters.  Further,  another  difficulty 
is  the  unavailability  of  a  high  fidelity  dynamic  model  for  hypersonic  vehicles.  This  is  largely 
because  of  the  inherent  difficulties  of  aerodynamic  modeling  in  hypersonic  flight  regime  as 
well  as  the  lack  of  the  sufficient  number  of  flight  tests  (limited  flight  tests,  that  have  been 
carried  out,  are  all  of  the  short  duration  and  within  a  very  restricted  flight  envelope  as  well) 
[Duan  and  Li,  2012,  Mirmirani  et  al.,  2005a].  These  uncertainties  in  aerodynamic  parameters 
require  a  control  system  design  based  on  robust  adaptive  control  techniques  that  learn  the 
uncertainty  on  the  data  and  provide  a  controller  that  is  able  to  maintain  desirable  flight 
performance  of  AHV 

Hence,  in  this  report,  two  kinds  of  robust  adaptive  control  design  approaches  are  pre¬ 
sented  for  the  stabilization  of  AHV  in  cruise  phase.  In  the  first  approach,  DI  technique  is 
used  to  derive  the  controller.  Further,  in  second  technique.  Barrier  Lyapunov  function  is  used 
to  derive  the  SC-GDI  technique.  In  AHV,  the  angle  of  attack  and  sideslip  angle  need  to  be 
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regulated  very  tightly  with  limited  control  surface  deflection  to  prevent  inlet  unstart  and  sub¬ 
sequently  vehicle  loss.  Hence,  both  control  techniques  are  augmented  with  newly  developed 
constrained  neuro- adaptive  technique  to  learn  the  unknown  disturbance  function  as  well  as 
its  Jacobian  matrix  very  quickly  Moreover,  BLF  is  used  to  formulate  constrained  weight  update 
law  for  unknown  disturbance  function  learning.  The  effectiveness  of  the  derived  controllers 
are  shown  by  promising  simulation  results. 


1.2  Literature  Survey 

Air-breathing  hypersonic  vehicle  (AHV)  research  gain  a  lot  of  attention  from  the  researchers 
because  of  its  explicit  advantages  over  conventional  air  vehicles.  Many  researchers  proposed 
various  models  of  AHV  and  shown  the  effectiveness  of  the  control  algorithm.  Since  the  AHV 
model  is  highly  uncertain  and  nonlinear,  hence,  robust  adaptive  nonlinear  controllers  are 
preferred  over  linear  controllers  for  controlling  of  AHV 

Longitudinal  Winged-cone  hypersonic  model  is  used  by  Marrison  and  Stengel  [1998]  to 
design  a  robust  control  by  assuming  parameter  uncertainties.  Marrison  and  Stengel  [1998] 
used  the  genetic  algorithm  to  find  unknown  parameter  values  by  selecting  stochastic  cost 
function.  Further,  a  combination  of  dynamic  inversion  and  the  stochastic  cost  function  is 
selected  with  the  genetic  algorithm  by  Wang  and  Stengel  [2000]  to  design  a  robust  controller. 
The  robust  adaptive  controller  is  proposed  by  [Xu  et  al.,  2012]  with  the  application  of  neural 
network  to  control  the  longitudinal  winged-cone  hypersonic  model. 

Another  longitudinal  hypersonic  vehicle  model  is  developed  by  Bolender  and  Doman 
[2005]  using  hypersonic  piston  theory.  Further,  various  practically  useful  variants  of  this  model 
are  proposed  in  the  literature.  Hypersonic  vehicle  dynamics  during  inlet  unstart  is  investigated 
by  Bolender  et  al.  [2009]  using  CFD.  Further,  nonlinear  dynamic  inversion  controller  design 
for  a  non-minimum  phase  hypersonic  vehicle  is  derived  by  Kuipers  et  al.  [2008].  Moreover, 
integrated  guidance  and  control  for  longitudinal  dynamics  of  the  vehicle  is  shown  by  Serrani 
et  al.  [2009] .  The  adaptive  controller  is  proposed  by  Gibson  and  Annaswamy  [2008]  in  presence 
of  thrust  and  actuator  uncertainties. 

The  third  type  of  hypersonic  longitudinal  model  is  developed  and  verified  by  Mirmi- 
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rani  et  al.  [2005a].  Further,  aero-propulsive  coupled  analysis  is  shown  by  Clark  et  al.  [2006]. 
Moreover,  Time-varying  vehicle  model  analysis  and  flexible  mode  analysis  of  the  vehicle  is 
presented  by  Fidan  et  al.  [2006].  Further,  Butt  et  al.  [2010],  Xu  and  Zhang  [2015]  formulated 
robust  adaptive  backstepping  control  methods  for  uncertain  AHV 

The  winged-cone  model  of  an  air-breathing  hypersonic  vehicle  is  developed  by  NASA  for 
the  analysis  of  six  degree-of-freedom  (six-DOF)  motion  of  equations.  Shaughnessy  et  al.  [1990] 
estimated  force  and  moment  coefficients  using  hypersonic  local  surface  inclination  code  and 
obtained  as  the  discrete  data  point.  Further,  these  aerodynamic  coefficients  are  plotted  as 
a  function  of  Mach  number,  the  angle  of  attack  and  control  surface  deflections.  Moreover, 
the  variation  of  Inertial  components  is  also  shown  wrt.  the  center  of  gravity  and  mass  of  the 
vehicle.  Further,  [Keshmiri  et  al.,  2005,  2004]  used  the  plots  of  discrete  data  points  (given  by 
Shaughnessy  et  al.  [1990])  and  obtained  the  expressions  of  aerodynamic  coefficients  using 
higher  order  interpolation.  Furthermore,  the  expressions  of  aero-coefficients  are  validated 
by  CFD  analysis  [Keshmiri  et  al.,  2006].  Six-DOF  motion  of  equations  for  AHV  are  derived  by 
Keshmiri  et  al.  [2007]  and  further  enhanced  to  ten  degree-of-freedom  motion  of  equations 
[Colgren  et  al.,  2009] .  An  adaptive  inverse  controller  for  AHV  is  proposed  by  Qian  and  liaxue 
[2010]  with  the  help  of  the  neural  network. 

Other  six-DOF  model  of  AHV  are  presented  by  Vick  [2014],  Chamitoff  [1992].  The  adap¬ 
tive  robust  controller  is  derived  by  Rollins  et  al.  [2013],  using  dynamic  inversion  technique  for 
the  model  proposed  by  Vick  [2014] .  Further,  the  same  model  is  used  by  Famularo  et  al.  [2016] 
to  obtain  a  state-constrained  control  solution  for  AHV  using  projection  operator.  Control- 
oriented  modeling  of  AHV  is  done  by  Huifeng  et  al.  [2011]  using  parameterized  configuration 
approach. 

feedback  linearization  concept  of  control  design  is  presented  by  Isidori  [2013]  for  various 
class  of  nonlinear  systems.  In  this  method,  the  nonlinear  system  is  transformed  in  a  way  such 
that  it  appeared  linear  in  transformed  coordinates.  Further,  using  stable  linear  error  dynamics, 
dynamic  inversion  control  technique  is  derived  by  Slotine  et  al.  [1991].  However,  this  control 
design  technique  assumed  that  the  system  is  perfectly  known,  which  is  very  rare.  Hence,  to 
overcome  this  limitation,  Padhi  et  al.  [2007]  presented  a  neuro- adaptive  dynamic  inversion 
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approach  for  robust  control  design,  in  which  the  bounded  disturbance  function  is  learned  by 
the  controller  over  time.  Further,  Ambati  and  Padhi  [2016]  augmented  neuro- adaptive  control 
design  technique  with  Jacobian  matrix  learning  of  unknown  disturbance  function  to  minimize 
the  transients. 

Furthermore,  Tee  et  al.  [2009]  formulated  the  output  constraint  nonlinear  control  for  strict- 
feedback  SISO  nonlinear  system  by  selecting  barrier  Lyapunov  function  (BLF)  as  a  Lyapunov 
candidate.  BLF  is  obtained  by  constraining  the  radially  unbounded  quadratic  Lyapunov 
function  (QLF)  into  a  locally  bounded  region.  Moreover,  backstepping  control  design  techni¬ 
que  used  by  Tee  et  al.  [2009]  to  obtain  output  constraint  control  for  both  known  and  uncertain 
SISO  nonlinear  systems.  Output  constrained  control  design  technique  is  further  enhanced  by 
selecting  time  varying  bound  for  uncertain  SISO  nonlinear  system  by  Tee  and  Ge  [2011] . 


1.3  Control  Objectives 

The  main  control  objective  in  this  report  is  to  design  a  robust  adaptive  control  algorithm 
that  stabilizes  the  AHV  during  the  cruise  flight  and  prevent  inlet  unstart.  Moreover,  the 
controller  should  be  robust  enough  to  take  care  of  model  uncertainties  and  disturbance, 
hence,  additional  controller  design  objective  is  to  minimize  the  need  of  high  fidelity  AHV 
model.  The  mission  scenario  is  shown  in  Fig.  1.1,  where,  the  AHV  is  required  to  follow  the 
desired  path  without  much  deviation.  Therefore,  further,  control  objectives  are  to  maintain 
the  states  as  well  as  the  error  in  disturbance  learning  within  the  prefixed  desired  bounds. 
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III  -  Cruise  Phase 


Figure  1.1  -  Mission  Scenario 
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Hypersonic  Vehicle  Model  Selection 
and  Vehicle  Dynamics 


2.1  Vehicle  Dynamics 

Rigid  winged-cone  model  of  air-breathing  hypersonic  vehicle  (AHV)  is  selected  in  this  report 
to  derive  the  robust  adaptive  controller.  The  AHV  geometry  is  shown  in  Fig.  2.1,  in  which 
propulsive  system  is  integrated  with  the  vehicle.  The  selected  vehicle  model  is  having  canard 
configuration,  however,  the  canard  is  used  only  for  subsonic/ supersonic  speed  regime.  As 
the  AHV  speed  approached  to  hypersonic  speed,  the  canard  fins  are  closed  inside  the  vehicle 
frame,  hence,  canards  have  no  effect  in  hypersonic  flight. 

The  vehicle  dimension  and  parameter  values,  obtained  by  Shaughnessy  et  al.  [1990], 
are  given  in  foot-pound-second  unit  system,  therefore,  first  values  are  converted  into  meter- 
kilogram- second  unit  system.  The  aerodynamic  coefficients,  presented  by  Keshmiri  et  al. 
[2005],  are  dimension  less  hence,  there  is  no  need  to  maintain  dimension  compatibility. 
Further,  the  important  vehicle  related  constants  used  in  controller  design  and  analysis  are 
summarized  in  Table  2.1 

The  vehicle  configuration  shown  in  Fig.  2.1  is  assumed  to  have  the  thrust  vector  orienta¬ 
tion  perfectly  aligned  with  body  x  axis.  Hence,  there  is  no  thrust  force  component  in  body  y 
and  z  directions. 


7 
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Figure  2.1  -  Air-breathing  Hypersonic  Vehicle  Geometry  [Shaughnessy  et  al.,  1990] 


Table  2.1  -  Important  AHV  Constants  Used  in  Controller  Design 


AHV  Constants 

Value 

Mass  (Mp) 

136080  kg 

Wing  Reference  Area  (S) 

334.7295 

Wing  Span  (h) 

18.2880  m 

Mean  Aerodynamics  Chord  (c) 

24.384  m 

Moment  Reference  Center  [Xmrc) 

37.7952  m 

Gravity  (g) 

9.8 

sec^ 

Maximum  Thrust  {Tmax) 

1467900  N 

Specific  Impulse  Dsp) 

1000  g  sec 

The  translational  dynamic  equations  of  AHV  are  given  as 


u  =  Rv-Qw-gsm6  + 


Pax  +  Ptx 


My 


V  =  wP  -  vR  + gcosO + 


^  ay 


(2.1) 


Faz 

w  =  uQ-  vP  +  gcosO  + - 


The  translational  dynamic  equations  (2.1)  of  AHV  can  also  be  transformed  into  Vf,  d, 
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and  j6  variables  as 


d 

$ 


UU+VV+  IVW 

uw-  wu 
+  iv'^ 

[vVt-vVt]  cos  j6 


(2.2) 


where,  total  AHV  velocity  (y^),  angle  of  attack  (a),  and  sideslip  angle  (j6)  is  calculated  as 

Vt  =  V  a=tan~^[^)y  p  = 

The  rotational  dynamic  equations  of  AHV  are  given  as 
P  =  CiQR  + CzPQ  +  C^Lam~^ 

Q  =  C5PR-C6{P^-R^)  +  C7Ma  (2.3) 

R  =  CgPQ  -  C2QR  + C4Lam~^ 

The  rotational  kinematic  equations  of  AHV  are  given  as 

(p  =  P  +  QsincptanO  +  Rcos(ptan9 

6  =  Qcoscp-  Rsincp  (2.4) 

y/  =  QsincpsecO  +  RcoscpsecO 

The  translational  kinematic  equations  of  AHV  are  given  as 

X  =  u  cos  6  cos  y/  +  (sin (p  sin 6  cos  t//^  -  cos  0  sin t//^)  i/  +  (cos  cp  sin 6  cos  t//^  +  sin 0  sin y/)  w 
y  =  ucos6smy/+  ism(psm6smy/  +  cos(pcosy/)v+  (cos 0 sin 0 sin t//^- sine/) cos t//^)u;  (2.5) 
h  =  usinO  -  vsincpcosO  -  wcoscpcosO 


Here,  w,  Vy  and  w  are  the  body  velocities  in  x,  y,  and  z  directions  respectively.  Further, 
roll,  pitch  and  yaw  angles  are  denoted  as  epy  6  and  y/.  P,  Q  and  R  are  roll  rate,  pitch  rate  and 
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yaw  rate  around  body  x,  y,  and  z  axis  respectively  Moreover,  h  is  the  altitude  of  the  vehicle  in 
inertial  frame  and  x  and  y  are  the  distances  covered  by  AHV  in  inertial  frame  of  reference. 

Furthermore,  Fax,  Fay,  and  Faz  are  the  aerodynamic  body  forces  acted  in  x,  y,  and  z 
directions  respectively.  Further,  Fj^  is  the  thrust  force  along  body  x  axis.  Moreover,  is  the 
mass  of  the  AHV  [Shaughnessy  et  al.,  1990]  and  g  is  the  gravity  force  acting  on  the  vehicle  in 
inertial  frame  of  reference.  The  moments,  Lam,  ^a,  and  Na  are  the  aerodynamic  moments 
around  center  of  gravity.  Further,  coefficients  Ci,  C2, . . . ,  Cg  are  calculated  from  inertia  values  of 
the  AHV  and  calculated  as 


Cl 

1 

1 

C2 

Ixzilzz  +  I XX  ~  lyy) 

C3 

1 

Izz 

C4 

[ixxhz  -  Ixz) 

lyz 

^8 

Ixx(Ixx-  lyy)  Ixz 

cg 

Ixx 

^5  —  I zz  Ixx) / lyy 
^6  “  ^xz!  lyy 

C'j  —  1/  lyy 

The  inertia  values  of  AHV  varies  with  respect  to  mass  because  mass  consumption  is 
assumed  during  flight.  Th  expression  for  the  inertia  values  are  obtained  by  Zhang  et  al.  [2010] 

-7.8809  X  10“®  Ml  +  25.88576  My  -  6.9683  x  10® 

-8.2890  X  10“^  Ml  +  265.9889  My  -  7.3048  x  10® 

lyy 

0 

where,  Ixx,Iyy,  and  Izz  are  diagonal  terms  of  inertia  matrix  and  Ixz  is  cross  diagonal 


and  given  as 

Ixx  — 

lyy  = 

Izz  — 
Ixz  — 
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term.  Further,  in  AHV;  the  fuel  consumption  is  very  high,  hence,  the  mass  variation  of  AHV  is 
considered  while  designing  the  controller.  The  mass  variation  dynamics  is  given  as 


My  =  (2.6) 

hpS 

where,  new  control  variable  apis  defined  as  aT=  Here,  T  is  the  required  thrust 

max 

magnitude  and  Tmax  is  the  maximum  thrust  that  AHV  can  generate  [Keshmiri  et  al.,  2004] .  The 
specific  impulse  of  thrust  is  denoted  by  Isp.  Further,  the  hypersonic  speed  generally  expressed 
in  therm  of  mach  number.  Hence,  the  conversion  from  velocity  to  mach  number  is  given  as 

Vt 

M  =  —  (2.7) 

a 

where  a  is  the  speed  of  sound  which  is  constant  for  cruise  flight  and  Vf  is  the  total 
velocity  of  the  vehicle. 


2.2  Forces  and  Moments  Calculation 

In  this  section,  the  calculation  of  forces  and  moments  used  in  six  degree-of-freedom  dynamics 
of  AHV  are  presented.  The  expressions  for  force  and  moment  coefficients  are  obtained  by 
Keshmiri  et  al.  [2005]  using  higher  order  interpolation.  Moreover,  the  obtained  aerodynamic 
coefficient  expressions  are  function  of  M,  a  (in  degree)  and  control  surface  deflections  (in 
degree).  Further,  the  AHV  consist  of  four  controls  5^,  5  a,  5r  and  ap-  Among  the  four  control 
variables,  5^,  5 a,  and  Sr  are  control  surface  deflection  for  right  elevon,  left  eleven  and  rudder 
respectively  and  ur  is  used  to  maintain  the  thrust  requirement  to  regulate  the  AHV  speed. 

The  calculation  of  aerodynamic  forces  are  done  in  the  similar  way  as  explained  in 
[Shaughnessy  et  al.,  1990]  to  maintain  the  compatibility.  The  drag  force  D  is  calculated  as 

D  =  qSCo  (2.8) 
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where,  'q  is  the  dynamic  pressure  obtained  as  ^  Vf  is  total  AHV  velocity  and 

p  is  the  pressure.  The  surface  area  of  the  vehicle  is  denoted  as  S.  The  drag  coefficient  Cols 
calculated  as 


Cd  -  Coa  +  CD,6e  ^D,5a  +  ^D,5r 


here,  Cdu  is  the  basic  drag  incremental  coefficient  and  ^D,6a  ^D,6r  ^^e  drag 

incremental  coefficients  for  right  elevon,  left  eleven  and  rudder  respectively 

Further,  the  Lift  force  L  is  calculated  as 


L  =  qSCi 


(2.9) 


where,  the  lift  coefficient  Cl  is  obtained  as 


Cl  -  Clu  +  CL,de  + 


here,  Clu  is  the  basic  lift  incremental  coefficient  and  Cl, 5^  and  Cl, 5^  are  lift  incre¬ 

mental  coefficients  for  right  elevon,  left  elevon  and  rudder  respectively. 

The  side  force  Y  is  calculated  as 


F  =  qSCy 


(2.10) 


where,  Cy  is  the  side  force  coefficient,  calculated  as 


Cy  -  CvfiP  +  CY,6e  +  ^Y,6a  +  CY,6r 


here,  Cy^  is  the  basic  vehicle  sideslip  derivative  and  j6  is  known  as  sideslip  angle  which 
is  in  radians.  Further,  CY,Sa  incremental  coefficients  for  right  elevon,  left 

elevon  and  rudder  respectively. 
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The  rolling  moment  L^c  is  given  as 

LAc  =  qSbCi^  (2.11) 

here,  C/„  is  the  rolling  moment  coefficient,  calculated  as 

(Pb\  (Rb\ 

Clm  -  Clpfi  +  Q,5,  +  Cl, 5a  +  Cl, 5,  +  Ci^  1^  — j  +  Ci^  1^  — j 

Here,  C/^  is  the  basic  rolling  moment  with  sideslip  derivative  and  are 

rolling  moment  incremental  coefficients  for  right  elevon,  left  eleven  and  rudder  respectively. 
Further,  C/^  and  C/^  are  rolling  moment  derivative  for  roll  rate  and  yaw  rate. 

The  pitching  moment  Mac  relative  to  moment  reference  center  is  calculated  as 

MAc  =  qScCm  (2.12) 


where.  Cm  is  the  pitching  moment  coefficient,  obtained  as 


Cm  —  Cma  +  Cm,6e  Cm, 5 a  Cm,6r  Cn 


Here,  Cma  is  the  basic  pitching  moment  incremental  coefficient  relative  to  moment 
reference  center  and  Cm,6,y  Cm,5ay  Cm,5r  ^re  pitching  moment  incremental  coefficients  for  right 
elevon,  left  elevon  and  rudder  respectively.  Moreover,  Cm^  is  pitch  rate  dynamic  derivative. 

The  yawing  moment  Nac  relative  to  moment  reference  center  is  calculated  as 

NAc  =  qSbCn  (2.13) 


where,  yawing  moment  coefficient  Cn  is  calculated  as 
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Here,  Cn^  is  the  basic  yawing  moment  with  sideslip  derivative  and  are 

yawing  moment  incremental  coefficients  for  right  elevon,  left  eleven  and  rudder  respectively. 
Further,  and  Cn,  are  yawing  moment  derivative  for  roll  rate  and  yaw  rate. 

The  aerodynamic  and  thrust  forces  used  in  six-DoF  equations  (2.1)  are  obtained  as 

Fax  =  Lsina-Dcosa 
Fay  =  y 

Faz  =  -Dsina-Lcosa 

Ft  X  —  ^  Tmax 

The  aerodynamic  moments  used  in  six-DoF  equations  (2.3)  are  calculated  as 

Fam  —  Fj{c 

—  ^Ac~^cgFaz 
Ffa  —  F[Ac~^^cgFay 


here,  M^,  and  Na  are  the  aerodynamic  moments  around  center  of  gravity.  Moreo¬ 
ver,  Lao  Macj  and  Nac  are  the  aerodynamic  moments  around  moment  reference  center  of 
vehicle.  Further,  theh  mass  of  the  AHV  is  veries  with  time,  hence,  center  of  gravity  Xcg  is  also 
varies  with  time.  The  expression  of  Xcg  (in  meter )  wrt.  mass  is  given  as  (Zhang  et  al.  [2010]) 

Xcg  =  1.6210  X  10“^'’ My -5.5226  X  10“^  M„  + 7.3578  (2.14) 
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Nominal  Dynamic  Inversion  Control 
Design  of  Hypersonic  Vehicle 


In  this  section,  three-loop-control  design  architecture  is  presented  for  the  air-breathing  hyper¬ 
sonic  vehicle  (AHV)  as  shown  in  Fig.  3.1.  Dynamic  inversion  technique  is  used  to  obtain  the 
controller  for  tracking  the  reference  signal.  Three-loop-control  design  consist  of  three  loops, 
the  outer  loop  is  used  to  obtain  the  desired  pitch  angle  0* .  Further,  the  desired  pitch  angle  0* 
is  the  angle  required  to  track  the  constant  altitude.  Furthermore,  the  output  of  the  outer  loop 
is  used  as  the  reference  command  for  the  intermediate  loop.  Hence,  by  using  DI,  the  desired 
roll  rate  P* ,  pitch  rate  Q* ,  and  yaw  rate  P*  are  obtained  to  track  the  reference  input.  Finally,  in 
the  inner  loop,  the  control  surfaces  deflection  is  calculated  to  track  the  desired  angular  rates 
obtained  from  the  intermediate  loop.  Further,  the  velocity  controller  is  used  to  maintain  the 
desired  velocity  of  the  vehicle  to  prevent  the  inlet  unstart  (low  vehicle  speed  is  also  a  reason  of 
inlet  unstart).  Further,  velocity  controller  is  used  to  obtain  the  required  amount  of  thrust,  that 
is  needed  to  maintain  the  desired  velocity  of  the  AHV  The  dynamic  inversion  (DI)  formulation 
of  all  three  loops  is  given  in  further  sections 


3. 1  Outer  Loop  Control  design 

The  controller  is  designed  for  cruise  phase  of  the  flight,  hence,  the  constant  altitude  is  reference 
signal  to  the  outer  loop.  Dynamic  inversion  is  used  to  obtained  the  desired  pitch  angle  0* 
from  h  equation  (2.5).  Further,  the  the  derivative  of  the  desired  reference  signal  is  assumed  to 


15 
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Figure  3.1  -  Block  Diagram  of  Three-Loop-Control  Design 


be  zero.  Hence,  using  dynamic  inversion,  the  error  dynamics  is  written  as 


eh  +  kheh  =  0 


here  eh  =  -  his  the  error  in  reference  height  and  obtained  height  and  kh  is  positive  constant 

gain.  Further,  the  error  dynamics  is  expended  as 

-h)  +  kh{h^ -h)  =  0 
h  =  kh[h^-h) 


Using  h  equation  from  equation  2.5,  the  error  dynamics  can  be  rewritten  as 
usinO  -  {vsmcf)+  wcosc[))  cosO  =  khih"^  -  h) 


On  solving  for  pitch  angle  (0  =  0*)  and  rearranging  the  terms,  the  desired  pitch  angle  0*  is 
obtained  as 


0*  =  tan 


( vsmcf)+ wcoscf)\  . 

khih*-h) 

1  u  I 

^\l  [vsin(p+  wcoscpY'  ^ 

(3.1) 


Furthermore,  desired  roll  angle  0*  and  desired  yaw  angle  t/A*  is  assumed  to  be  zero  as  reference 
values  for  intermediate  loop. 
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3.2  Intermediate  Loop  Control  design 

The  output  of  the  outer  loop  is  act  as  the  reference  signals  for  the  intermediate  loop.  The 
desired  roll  rate  P*,  pitch  rate  Q*,  and  yaw  rate  P*  are  obtained  by  using  Dl.  The  error 
dynamics  of  the  intermediate  loop  is  written  as 


o 

o 

ee 

+ 

o 

o 

ee 

^yr 

0  0 

eyr 

Here,  =  0*  -  0,  is  the  error  between  the  desired  roll  angle  and  obtained  roll  angle.  Similarly, 
=  0*  -  0,  is  the  error  between  desired  and  obtained  pitch  angle  and  ey/  =  -  xf/,  is  the 

error  between  desired  and  obtained  yaw  angle.  Moreover,  A:^,  Icq  and  ky;  are  positive  definite 
constants.  Further,  rearranging  the  intermediate  loop  error  dynamics  using  equations  2.4,  the 
desired  body  rates  are  obtained  as 


p* 

1  tan  0  sin  0  tan  0  cos  0 

-1 

kcpifp’" -(p) 

Q* 

0  cosc[)  -sincf) 

keie*-e) 

P* 

0  sin(/)sec0  coscf)sec6 

k^,  {y/*  -  y/) 

3.3  Inner  Loop  Control  design 

In  this  section,  the  control  surfaces  deflection  are  obtained  as  an  output  using  Dl  on  the  inner 
loop.  In  the  inner  loop,  the  desired  body  rates  obtained  from  intermediate  loop  are  considered 
as  reference  signals  and  equations  2.3  are  used  to  formulate  the  error  dynamics.  Further,  the 
inner  loop  error  dynamics  is  written  as 


ep 

o 

o 

ep 

eQ 

+ 

0  kq  0 

eQ 

ep 

o 

o 

ep 
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Here,  ep  =  -  P,  is  the  error  between  the  desired  roll  rate  and  obtained  roll  rate.  Similarly, 

eg  =  Q*  -  Q,  is  the  error  between  desired  and  obtained  pitch  rate  and  Cr  =  -  R  is  the 

error  between  desired  and  obtained  yaw  rate.  Moreover,  kp,  kq  and  kp  are  positive  definite 
constants.  Further,  rearranging  the  error  dynamics  using  equations  2.3,  the  following  nonlinear 
coupled  equations  are  obtained  as 


/Cp(P-P*)  +  C3Lam=0 
C5PP+/Cq(Q-Q*)  +  C7M^=0 
c^PQ  +  kR[R-R^)  +  c^Na=^ 


The  above  equations  contained  only  control  surface  deflections  5^,  5 ay  and  5r  as  an  unknown. 
Hence,  Newton-Raphson  algorithm  is  used  to  obtain  the  desired  control  surface  deflections. 

3.4  Velocity  Loop  Control  design 

Velocity  control  loop  is  an  extra  loop  presented  in  the  control  design  to  maintain  the  constant 
desired  velocity  of  the  AHV  The  forward  velocity  u  is  the  major  component  of  the  AHV  velocity 
and  other  velocity  components  are  negligible  as  compared  to  forward  velocity.  Hence,  dynamic 
inversion  technique  is  used  on  forward  velocity  equation  given  in  equation  2.1  to  obtained 
the  desired  thrust  magnitude  for  maintaining  the  desired  velocity.  The  velocity  loop  error 
dynamics  is  written  as 


—  0 


here  Cu  =  -  u,  is  the  error  between  the  desired  froward  velocity  and  obtained  forward 

velocity  and  ku  is  positive  constant  gain.  Further,  the  error  dynamics  is  expended  as 


ii  =  kuiu"  -  u) 


Further,  the  desired  ap  is  obtained  using  equation  2.1  as 


(3.2) 
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3.5  Simulation  Results 

3.5. 1  Calculation  of  Trim  Conditions 

Trim  conditions  needed  to  be  calculated  before  carrying  out  the  simulations.  Trim  con¬ 
ditions  are  the  the  conditions  around  which  the  aircraft  is  in  equilibrium.  To  obtain  trim 
conditions  for  the  cruise  flight,  the  constant  height  and  constant  forward  velocity  are  as¬ 
sumed  as  /z*  =  24000  m  and  w*  =  1786.4  m/s.  Accordingly,  P,  Q,  R  and  (p  is  assumed  as 
zero.  Further,  the  variables  [iz,  w,  0,  5^,  6 ay  Sry  (Tt]^  can  be  calculated  from  the  equations 
[iiy  [>,  iF,  P,  Q,  R,  h]^  =  0  using  Newton-Raphson  method.  Further,  t/z,  x  and  y  are  not  ap¬ 
peared  in  any  of  the  equations  used  for  trim  calculation,  hence  y/,  x  and  y  equations  are 
not  included  in  trim  calculation.  Finally,  the  trim  conditions  are  obtained  as  iz  =  0  m/5, 
w  =  91.8  m/5,  6  =  2.94°,  Se  =  -7.21°,  6a  =  -7.21°,  5,  =  0°,  or  =  0.194. 


3.5.2  Simulation  results 

In  this  section,  the  results  obtained  from  the  nominal  controller  is  presented.  The  results  are 
obtained  by  varying  the  initial  conditions  around  the  trim  values  to  show  the  applicability  of 
the  obtained  controller.  For  simulations,  rk^  method  is  used  with  step  size  as  0.01 5cc.  The 
results  are  simulated  for  150  sec.  All  the  vehicle  parameter  and  force-moment  coefficients  are 
taken  from  the  literature  published  by  [Shaughnessy  et  al.,  1990]  and  Keshmiri  et  al.  [2005]. 

The  initial  conditions  around  the  trim  values  are  given  as,  u  =  1836.4  m/5,  v  =  0  m/s, 
w  =  91.8  m/5,  P  =  3  deg/ sec,  Q  =  2  deg/ sec,  R  =  -3  deg/ sec,  (p  =  0  deg,  6  =  2.94  deg, 
y/  =  0  deg,  x  =  Qm,y  =  Qm,h  =  24100  m  and  initial  mass  =  136077.6  kg.  The  inner  loop 
gains  kp,  kq,  and  kp  are  selected  as  2.0  and  the  intermediate  loop  gains  /c^,  kg,  and  ky;  are 
selected  as  0.2.  Further  the  outer  loop  control  gain  is  selected  as  0.01  and  the  velocity  loop 
controller  gain  is  chosen  as  0.01. 

The  velocity,  angle  of  attack  and  sideslip  angle  profiles  are  shown  in  Fig.  3.2  respectively. 
Velocity  is  perfectly  tracked  by  the  application  of  derived  nominal  controller.  The  angle  of 
attack  is  within  the  bound  to  prevent  the  inlet  unstart.  Further,  sideslip  angle  is  zero  after 
sometimes  hence  the  reference  trajectory  is  effectively  tracked  by  the  AHV 
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Figure  3.2  -  Mach  Number  (M),  Angle  of  Attack  (a)  and  Sideslip  Angle  (j6)  Profiles  of  AHV 
using  Nominal  Dynamic  Inversion 


Further,  the  profiles  of  components  of  the  velocity  are  shown  in  Fig.  3.3  respectively.  The 
forward  velocity  u  is  the  major  component  of  velocity  and  become  constant  after  a  finite  time, 
hence  the  required  velocity  is  maintained  by  the  AHV  to  prevent  inlet  unstart.  Further,  the 
cross  velocity  u  become  zero,  hence  there  is  side  deviation.  Furthermore,  the  altitude  profile 
is  shown  in  Fig.  3.4.  The  reference  altitude  is  nearly  tracked  by  the  states,  however,  tracking 
altitude  is  not  a  major  concern  here,  hence,  the  results  are  acceptable.  The  continuous  and 
smooth  mass  and  flight  path  angle  profiles  of  the  AVH  are  shown  in  Fig.  3.4.  The  flight  path 
angle  is  become  zero  after  finite  time,  which  means  the  AHV  is  flying  at  zero  angle  wrt.  inertial 
frame  of  reference. 

Roll,  pitch,  and  yaw  angle  profiles  are  presented  in  Fig.  3.5.  The  roll  and  yaw  angle 
profiles  are  perfectly  tracking  zero  in  spite  of  the  variation  in  initial  conditions  hence,  the  AHV 
is  stabilized  after  the  control  application  (AHV  is  inherently  unstable).  Further,  the  pitch  angle 
profile  is  stabilized  at  a  non-zero  value,  which  is  essential  to  generate  the  necessary  lift  in  the 
AHV 
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Figure  3.3  -  Velocity  Profiles  w,  Vy  w  in  x,  y,  and  z  Directions  Respectively  using  Nominal 
Dynamic  Inversion 


The  angular  rate  profiles  are  shown  in  Fig.  3.6  respectively.  After  some  transients,  the 
angular  rate  profiles  are  stabilized  at  zero,  which  means  AHV  reached  the  equilibrium  point. 
Hence,  the  controller  formulated  using  nominal  DI  is  effective  for  control  and  stabilization  of 
AHV  Furthermore,  the  stabilization  time  for  inner  loop  variables  is  lesser  than  the  intermediate 
loop  variables  because  of  the  three-loop-control  design  methodology.  The  control  profiles  are 
shown  in  Fig.  3.7.  The  control  surface  deflections  are  within  the  saturation  limits  [-40°,  30°] 
hence,  the  solution  provided  by  the  controller  is  physically  realizable.  Further,  the  control 
(Jr  is  saturated  within  the  bound  [0.1, 0.9],  hence,  the  velocity  profile  takes  times  to  track  the 
desired  trajectory.  However,  the  stability  of  the  vehicle  is  ensured  by  the  derived  nominal 
controller. 
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Figure  3.4  -  Altitude  (/z),  Mass  {My)  and  Flight  Path  Angle  (y)  Profiles  of  AHV  using  Nominal 
Dynamic  Inversion 
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Figure  3.5  -  Roll  Angle  (0),  Pitch  Angle  (0),  and  Yaw  Angle  (t/z)  Profiles  of  AHV  using  Nominal 
Dynamic  Inversion 
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Figure  3.6  -  Roll  Rate  (P),  Pitch  Rate  (Q),  and  Yaw  Rate  (R)  Profiles  of  AHV  using  Nominal 
Dynamic  Inversion 


Figure  3.7  -  Control  Profiles  using  Nominal  Dynamic  Inversion 
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Further,  the  simulation  results  are  presented  for  various  different  initial  conditions 
around  trim  values  to  show  the  effectiveness  of  the  nominal  D1  controller.  The  random 
perturbation  in  forward  velocity  u  is  assumed  in  between  [-100,  100]  mis  around  trim  value. 
Moreover,  the  initial  condition  random  perturbation  in  angular  rates  are  assumed  in  between 
[-5,  5]  deg!  sec  around  trim  values.  Further,  the  random  perturbation  in  altitude  is  assumed 
in  between  [-500,  500]  m  around  trim  value. 

Velocity,  angle  of  attack  and  sideslip  angle  profiles  are  shown  in  Fig.  3.8  and  Fig.  3.9  for 
various  initial  conditions.  The  velocity  profiles  for  each  initial  condition  follow  the  same  trend 
and  become  constant  after  some  time.  Similarly,  the  angle  of  attack  profiles  are  also  within  the 
accepted  range  to  prevent  inlet  unstart.  Further,  the  altitude  profiles  are  shown  in  Fig.  3.10  for 
various  initial  conditions.  All  the  altitude  profiles  are  converged  around  the  reference  value. 
Furthermore,  mass  and  flight  path  angle  profiles  are  shown  in  Fig.  3.10  for  perturbed  initial 
conditions.  The  flight  path  angle  profiles  are  converged  to  zero  within  the  finite  time,  which 
shown  the  effectiveness  of  the  derived  controller. 
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Figure  3.8  -  Mach  Number  (AT),  Angle  of  Attack  (a)  and  Sideslip  Angle  (j6)  Profiles  for  Different 
Initial  Conditions  using  Nominal  Dynamic  Inversion 
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Figure  3.9  -  Velocity  Profiles  w,  Vy  w  in  x,  y,  and  z  Directions  Respectively  for  Different  Initial 
Conditions  using  Nominal  Dynamic  Inversion 
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Figure  3.10  -  Altitude  (/z),  Mass  (My)  and  Flight  Path  Angle  (y)  Profiles  for  Different  Initial 
Conditions  using  Nominal  Dynamic  Inversion 
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Figure  3.11  -  Roll  Angle  (0),  Pitch  Angle  (0),  and  Yaw  Angle  (if/)  Profiles  for  Different  Initial 
Conditions  using  Nominal  Dynamic  Inversion 
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Figure  3.12  -  Roll  Rate  (P),  Pitch  Rate  (Q),  and  Yaw  Rate  (P)  Profiles  for  Different  Initial 
Conditions  using  Nominal  Dynamic  Inversion 
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Figure  3.13  -  Control  Profiles  for  Different  Initial  Conditions  using  Nominal  Dynamic  Inversion 


The  angular  angle  and  rate  profiles  are  shown  in  Fig.  3.11  and  Fig.  3.12  respectively  for 
various  initial  conditions.  All  the  state  trajectories  are  converged  to  desired  trajectories.  Hence, 
the  performance  of  designed  nominal  controller  remains  good  for  perturbed  initial  conditions 
and  the  AHV  is  stabilized  for  each  initial  condition.  Finally,  the  control  profiles  are  shown  in 
Fig.  3.13,  where  all  the  control  trajectories  are  smooth  and  within  the  saturation  bound. 
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Neuro- adaptive  Dynamic  Inversion 
with  Jacobian  Learning 


In  this  section,  the  formulation  of  adaptive  controller  using  neuro -adaptive  dynamic  inversion 
technique  is  presented  with  necessary  mathematical  details.  In  neuro -adaptive  DI  technique, 
first,  the  nominal  system  model  is  assumed  and  the  desired  controller  is  obtained  by  using  DI 
technique.  Further,  the  actual  model  is  approximated  into  an  approximated  system  model  and 
then  using  DI  technique,  the  adaptive  controller  is  designed  by  enforcing  the  error  between 
the  approximated  and  desired  model  to  zero.  The  approximation  of  the  actual  system  model 
and  formulation  of  adaptive  controller  is  explained  in  further  sections. 


4. 1  Neuro-adaptive  Dynamic  Inversion  Control  Formulation 

Consider  a  nonlinear  nominal  plant  with  completely  known  dynamics  as 

Xd  =  nXd.Ud)  {m<n)  (4.1) 

here,  is  the  nominal  control  obtained  from  nominal  dynamic  inversion  technique 
(chapter  3).  However,  in  presence  of  external  disturbance  or  parameter  uncertainty,  the 
nominal  controller  may  not  perform  satisfactorily,  hence,  the  adaptive  controller  is  needed  to 
maintain  the  performance  of  the  vehicle. 

The  adaptive  controller  is  formulated  for  the  actual  system  model  having  model  uncer- 


28 
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tainties  and  disturbances.  Hence,  consider  the  nonlinear  dynamics  of  the  system  as 

X  =  f{X,U)  +  d{X)  XeSR"",  {m<n)  (4.2) 

where  d[X)  is  the  external  disturbance  or  parameter  uncertainties  present  in  the  system. 
In  the  neuro- adaptive  technique,  first  the  actual  model  is  approximated  by  learning  the 
disturbances  presented  in  the  actual  model.  Further,  the  Lyapunov  method  is  used  to  obtained 
the  disturbance  learning  rule.  The  approximated  model  of  the  actual  system  is  considered  as 

=  /(X,  U)  +  d[X)  +  kap[X  -  XJ  X^  e  (4.3) 

here,  d{X)  is  the  approximation  of  the  disturbance  d(X),  present  in  the  actual  system. 
kap  is  the  positive  definite  matrix.  Further,  the  error  between  the  desired  system  and  approxi¬ 
mated  system  is  enforced  to  zero  by  using  linear  error  dynamics,  as 

Xa-Xd  +  kadiXa-Xd)  =  0  (4.4) 


here,  kad  is  the  positive  definite  matrix.  Further,  equation  4.4  is  rearranged  in  a  general 
form  as 


/(X,  X^,  Xd.  H,  Ud.  d{X))  =  0  (4.5) 

The  expression,  obtained  in  equation  4.5,  is  the  general  expression  to  obtain  the  control. 
In  the  expression,  the  adaptive  control  U  is  obtained  by  using  Newton-Raphson  method, 
provided  d{X)  is  known.  Further,  d{X)  is  obtained  by  weight  update  rule  explained  in  section 

4.2. 


4.2  Constrained  Weight  Update  Rule  with  Jacohian  Learning 

In  this  section  the  learning  rule  for  d(X)  is  derived,  where,  d(X)  is  the  approximation  of  the 
disturbance  d{X).  Further,  the  error  between  the  actual  and  approximated  system  is  expressed 
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as 


E  =  X-Xa 


(4.6) 


by  differentiating  equation  4.6,  following  is  obtained 


E  =  X-Xa 


(4.7) 


further,  substituting  equations  4.2  and  4.3  into  4.7,  the  resulting  expression  is  written  as 


E  =  diX)-diX)-kap{X-Xa) 


(4.8) 


further,  for  analysis  purpose,  the  error  is  decomposed  into  individual  channels  as  c;  = 
Xi  -  Xi^,  where,  i  =  l,2 . n.  Hence,  the  channel  error  dynamics  can  be  written  as 


ei  =  di[X)-di{X)-ki^^{ei) 


(4.9) 


The  disturbance  is  assumed  to  be  linear  in  parameter  form,  hence 


diiX)  =  wf(piiX)+€i 

diiX)  =  wJcPilX) 


here,  (pi  (X)  is  the  basis  function  and  Wi  is  the  ideal  unknown  weights.  Moreover,  Wi  is 
the  approximated  weights  and  c/  is  the  approximation  error. 

Further,  the  error  dynamics  (4.9)  can  be  written  as 


Ci  =  iuf  (pi  (X)  +  c,  -  (Ci) 


(4.10) 


where,  =  wi  -  Wj,  and  further,  =  -Wi.  Moreover,  to  obtain  the  learning  rule,  the 
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combination  of  BLF  and  QLF  is  selected  for  each  channel  as 


wT  Wi 
2ri 


+ 


ddiiX) 

ddiiX)' 

T 

ddiiX) 

ddiiX) 

dX 

dX 

2 

dX 

dX 

(4.11) 


here,  /c^.,  y/,  and  0/  are  the  positive  definite  quantities  and  ji  is  a  positive  constant 
known  as  learning  rate.  Moreover,  and  are  the  Jacobian  matrices  of  the  distur¬ 

bance  function  and  its  approximation,  which  are  used  to  minimize  the  transients.  Further,  the 
BLF  is  used  to  constrained  the  error  e/  within  a  constant  positive  bound  |/c^J. 

further,  taking  the  derivative  of  the  Lyapunov  function  (4.11)  and  rearranging  the  com¬ 
ponents,  following  is  obtained 


~  T 

w. 


ei(piiX)  Wi 

d(pi 

d(pi 

T 

kl-e]  Ji 

dX  _ 

dX  _ 

eiCi 


'Ci  e. 

^ap  I 


bi  i  hi 


,2  + 


d(pi 

d(pi 

dX  _ 

^'dt 

dX  _ 

Wi 


(4.12) 


here,  ^  is  the  Jacobian  matrices  of  the  basis  function.  Further,  by  equating  wf  =  0,  the 
weight  update  rule  is  obtained  as 


Wi 


Ip  ^  \d<pi 

Ji  \  dX 


0/ 


d(pi 

dX 


T 


^  ei(l>iiX) 


(4.13) 


here.  Ip  is  the  identity  matrix  of  dimension  p.  Further,  the  diiX)  is  obtained  by  the 
expression  di  (X)  =  wj (pi  (X) .  Moreover,  the  further  analysis  of  Vi  equation  (4.12)  can  be  found 
in  the  literature  published  by  Ambati  and  Padhi  [2016],  in  which,  the  boundedness  of  the  error 
is  shown  mathematically. 


4.3  Simulation  Results 

In  this  section,  the  simulation  results  are  presented  with  the  application  of  adaptive  controller 
on  the  actual  system  obtained  by  neuro-adaptive  DI  technique.  The  model  uncertainties  are 
assumed  in  the  aerodynamic  coefficients  and  in  the  inertia  values.  The  disturbance  is  added 
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in  P,  Q,  and  R  state  equations.  The  initial  conditions  are  given  as,  Vf  =  1789  m/5,  a  =  2.94  degy 
j6  =  0  degy  P  =  0  degIseCy  Q  =  0  degIseCy  R  =  0  degIseCy  (p  =  0  degy  6  =  2.94  degy  y/  =  0  degy 
x  =  t)  rUy  y  =  0  myh  =  24000  m  and  initial  mass  My  =  136077.6  kg.  The  inner  loop  gains  /cp, 
kQy  and  Icr  are  selected  as  2.0  and  the  intermediate  loop  gains  /c^,  /c^,  and  ky/  are  selected  as 
0.1.  The  outer  loop  control  gain  k^  is  selected  as  0.01  and  further,  the  velocity  loop  controller 
gain  is  chosen  as  0.1. 


The  disturbance  learning  rate  7/  is  selected  as  40  and  the  constraint  on  the  error  (equa¬ 
tion  4.6)  is  selected  as,  0.7  deg!  sec.  The  positive  definite  matrix  kad  used  in  equation  4.4  is 
selected  as  7,  5)  and  used  in  equation  4.3  is  selected  as  d/ag(15,  18,  15).  Here,  the 

gain  values  of  kap  is  more  then  kad  so  that  Xa  (very  quickly)  before  Xa^  Xd-  Further,  the  si¬ 
mulation  is  carried  out  far  random  disturbance  within  30%  of  model  uncertainty.  Furthermore, 
the  basis  functions  (/)(x)  is  selected  as  [l,  P,  P^,  Q,  Q^,  (PQ)^,  P,  P^,  (PP)^,  (QP)^,  (PQP)^] 


and  the  lacobian  matrix  of  the  basis  function 


is  calculated  as 


d(p 

dx 


0  1  2P  0  0  2PQ2  0  0  2PP2  0  2PQ2p2 

0  0  0  1  2Q  2P^Q  0  0  0  2QR^  2P^QR^ 

00000  0  1  2P  2P^R  2Q^R  2P^Q^R 


(4.14) 


Further,  the  velocity,  angle  of  attack  and  sideslip  angle  profiles  are  shown  in  Fig.  4.1 
respectively.  In  the  plots,  the  profiles  obtained  using  nominal  control  for  actual  system  are 
become  unbounded,  hence,  the  nominal  controller  is  not  effective  in  presence  of  uncertain¬ 
ties.  Further,  the  desired  trajectories  are  effectively  tracked  by  the  obtained  profiles  with  the 
application  of  adaptive  controller.  Moreover,  the  angle  of  attack  is  within  the  bound  and  the 
constant  velocity  is  maintained  by  the  velocity  controller  to  prevent  the  inlet  unstart. 

Similarly,  the  constant  height  is  maintained  by  the  adaptive  controller  as  shown  in 
the  Fig.  4.2.  The  mass  profile  is  smooth  and  continuous  and  the  flight  path  angle  profile  is 
stabilized  at  zero.  Further,  the  roll,  pitch,  and  yaw  angle  profiles  are  presented  in  Fig.  4.3 
respectively.  The  states  are  bounded  and  tracking  the  desired  states.  Moreover,  the  pitch  angle 
is  non-zero  to  achieve  the  required  lift  in  the  vehicle. 
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Time  [sec] 

Figure  4.1  -  Mach  Number  (M),  Angle  of  Attack  (a)  and  Sideslip  Angle  (j6)  Profiles  of  AHV 
using  Neuro-adaptive  Dynamic  Inversion 


The  angular  rate  profiles  are  shown  in  Fig.  4.4  respectively.  After  some  transients,  the 
angular  rate  profiles  are  tracking  the  desired  states  perfectly,  which  means  AHV  is  stable. 
Hence,  the  formulated  neuro-adaptive  D1  controller  is  effective  to  control  the  uncertain 
AHV.  Further,  the  control  profiles  are  shown  in  Fig.  4.5.  The  control  surface  deflections  are 
continuous  (zoomed  in  box)  and  within  the  saturation  limits  hence,  the  solution  provided 
by  the  controller  is  physically  realizable.  Further,  the  neuro-adaptive  D1  control  profiles  are 
different  from  the  nominal  D1  control  profiles  because  of  the  uncertainty  presents  in  the 
system. 
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Figure  4.2  -  Altitude  (/z),  Mass  (Mp)  and  Flight  Path  Angle  (y)  Profiles  of  AHV  using  Neuro- 
adaptive  Dynamic  Inversion 
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adaptive  Dynamic  Inversion 
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Figure  4.5  -  Control  Profiles  obtained  using  Neuro-adaptive  Dynamic  Inversion 
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Further,  the  simulation  results  are  presented  for  random  disturbances  bounded  between 
[-30%,  30%] .  The  velocity  and  angle  of  attack  profiles  are  shown  in  Fig.  4.6  for  the  AHV.  The 
inlet  unstart  condition  is  successfully  prevented  by  the  proposed  controller  because  the  angle 
of  attack  is  perfectly  bounded  between  [-5°,  5°]  and  velocity  level  is  maintained  accordingly. 
The  altitude  and  mass  profiles  for  AHV  are  shown  in  Fig.  4.7.  The  smooth  and  continuous 
profiles  are  obtained  in  random  disturbances  shows  the  applicability  of  the  adaptive  controller. 
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Figure  4.6  -  Mach  Number  (AT),  Angle  of  Attack  (a)  and  Sideslip  Angle  (j6)  Profiles  for  Random 
Disturbances  using  Neuro- adaptive  Dynamic  Inversion 


The  angular  angles  rates  profiles  are  shown  in  Fig.  4.8  and  Fig.  4.9  respectively  for  the 
AHV.  All  the  state  trajectories  converged  to  desired  trajectories  within  finite  time  with  very 
minimal  transients.  Hence,  the  designed  adaptive  controller  stabilized  the  AHV  in  presence 
of  model  uncertainties.  Finally,  the  control  profiles  are  shown  in  Fig.  4.10,  where,  all  the 
trajectories  are  smooth  and  achievable  by  control  surfaces.  The  magnitude  of  the  required 
control  is  depend  on  the  level  of  model  uncertainty. 
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Figure  4.7  -  Altitude  (/z),  Mass  (M^;)  and  Flight  Path  Angle  (y)  Profiles  for  Random  Disturbances 
using  Neuro-adaptive  Dynamic  Inversion 
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Figure  4.8  -  Roll  Angle  (0),  Pitch  Angle  (0),  and  Yaw  Angle  iy/)  Profiles  for  Random  Disturbances 
using  Neuro-adaptive  Dynamic  Inversion 
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Figure  4.9  -  Roll  Rate  (P),  Pitch  Rate  (Q),  and  Yaw  Rate  (P)  Profiles  for  Random  Disturbances 
using  Neuro-adaptive  Dynamic  Inversion 
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Figure  4.10  -  Control  Profiles  obtained  for  Random  Disturbances  using  Neuro-adaptive  Dyna¬ 
mic  Inversion 
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Neuro-adaptive  Generalized  Dynamic 
Inversion  with  Jacobian  Learning 


5.1  Neuro-adaptive  State-constrained  Generalized  Dynamic  Inver¬ 
sion  Formulation 

In  this  section,  first  the  state-constrained  generalized  dynamic  inversion  (SC-GDI)  formulation 
is  presented  for  controlling  the  nonlinear  system  given  in  equation  4.2.  In  the  formulation, 
barrier  Lyapunov  function  (BLF)  is  selected  as  the  lyapunov  candidate  for  incorporating  the 
state  constraint  on  the  states.  Further,  it  has  been  shown  that  the  dynamic  inversion  technique 
is  a  subset  of  SC-GDI  and  by  selecting  a  tuning  variable,  the  SC-GDI  is  converted  into  DI. 
Further,  it  has  been  shown  that  the  obtained  controller  is  robust  enough  to  maintain  the 
constraint  on  states  in  presence  of  disturbances,  however,  to  achieve  perfect  tracking,  SC-GDI 
is  augmented  with  disturbance  learning  mechanism.  Since  the  disturbance  learning  rule 
is  already  explained  in  section  5.3,  hence,  in  this  chapter,  only  the  formulation  of  SC-GDI 
technique  and  the  derivation  of  neuro-adaptive  SC-GDI  control  is  presented. 

5.1.1  Formulation  of  State-Constrained  Generalized  Dynamic  Inversion 

State-constrained  generalized  dynamic  inversion  formulation  is  shown  in  this  section  for 
an  ideal  system  (no  uncertainty).  In  this  section,  therefore,  general  nonlinear  system,  given 
in  equation  4.2,  is  considered  by  assuming  diX)  =  0.  Further,  is  has  been  assumed  that  the 
dimension  of  control  {U  e  states  (X  e  and  output  (F  e  are  same  i.e.,  p  =  m  =  n. 

39 
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This  assumption  is  relevant  in  AHV  control  study  because  the  disturbance  is  assumed  in 
Py  Q,  and  R  state  equations  (2.3),  hence,  the  output  (P,  Q,  P),  states  (P,  Q,  R)  and  controls 
(5g,  5 ay  5r)  are  equal.  Moreover,  without  loss  of  generality,  the  problem  is  assumed  as  a 
regulation  problem,  which  is  also  relevant  too  in  AHV  control  because,  the  desired  P^,  and 
R^  are  always  zero. 

Further,  the  error  equation  is  defined  as 

E  =  Y  (5.1) 


for  the  purpose  of  analysis,  the  error  (P  e  K'^)  equation  (5.1)  is  decomposed  into  indivi¬ 
dual  channels.  Hence,  the  the  error  equation  is  given  as 

=  Yjy  7  =  l,2,...,m  (5.2) 


Ci  =  yiy  i  =  m  +  lym  +  2y...yn 


(5.3) 


The  combination  of  barrier  Lyapunov  function  (BLF)  and  quadratic  Lyapunov  function 
(QLF)  is  selected  as  the  Lyapunov  function  candidate.  The  BLF  is  chosen  for  the  states  having 
constraint  [Xj)  and  QLF  is  selected  to  maintain  the  boundedness  of  the  remaining  states  x/. 
Hence,  the  Lyapunov  function  is  written  as 


V 


^Liog 


A;?,  -  e^; 


4  2 

i=m+l 


el 


(5.4) 


here,  ksj  is  the  positive  constant,  selected  as  the  constraint  on  ej  i.e.,  ej  <  \  ksj  |  V  r  >  0. 
Further,  the  derivative  of  the  Lyapunov  function  (5.4)  can  be  obtained  as 


V 


L 


ejej 


t.2  L  eiSi 

=  1  •^Sj  i=m+l 


(5.5) 
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further,  for  compact  analysis,  equation  5.5  is  written  into  vector  format  as 

V  =  [BiEfE  (5.6) 


where  E=lei,  62,...,  e„]^  and  [Bl]  e  is  a  symmetric  diagonal  matrix,  given  as 


[Bl]  =  diag\ 


*^1  '^$2  *^2 


ki  -  ei 


(5.7) 


By  using  Lyapunov  stability  theory,  for  a  stable  close  loop  system  (V  <  o).  The  stability 
of  the  system  is  enforced  by  equating  V  with  a  negative  definite  function,  defined  as 


V  =  -E^KbifE 


(5.8) 


where,  [Kbi /]  e  is  a  positive  definite  symmetric  diagonal  matrix  expressed  as 

ki  kfn 


[Kbif]=diag\ 


{ki-elf[kl-el,f 


y  km+ly  •  •  •  ykn 


here,  [iQ,/ f]  >0,  Vf  >  0  because  /c^.  >  e/,  V r  >  0  (Tee  et  al.  [2009])  and  /q,  /C2, . . . ,  kn  are 
selected  as  positive  constants.  Further,  A  e  5R  and  A  can  be  different  for  different  terms  in 
[Ki)i  f] ,  however,  for  simplicity  unique  value  of  A  has  been  selected  in  this  report.  Further,  by 
equating  the  equation  5.6  with  5.8,  one  can  obtain 


iBLEfE  =  -E^Ki,ifE 


(5.9) 


Since  bJ  =  Bi  hence,  equation  5.9  can  be  rewritten  as 

E^BLE  =  -E^KbifE 


(5.10) 
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Finally,  the  error  dynamics  can  be  obtained  as 


E  +  Bl^KhifE  =  0 


(5.11) 


where  [B^^K^if]  e  given  as, 


Here,  equation  5.11  is  known  as  the  nonlinear  error  dynamic  equation,  used  for  control 
calculation.  Further  for  control  calculation,  equation  5.11  is  rewritten  by  using  equation  5.1 
and  expressed  as. 


X  +  Bl^KhifX  =  0 


(5.12) 


Hence,  the  general  expression  for  obtaining  the  control  is  given  as 


f{X,U)+Bl^KbifX  =  0 


Hence,  by  using  the  Newton-Raphson  method,  the  control  U  is  obtained.  Further,  if  the 
system  dynamics  is  in  control  affine  form  i.e.,  X  =  /(X)  +  g(X)  U,  where,  /(X)  and  g(X)  are 
nonlinear,  smooth  function  of  X,  then  the  state-constrained  control  solution  U  is  obtained  as 


[g{X)]-^{-Bl^KMfX-f{X)) 


u 


(5.13) 


Note.  By  careful  observation,  it  can  be  seen  that  by  selecting  A  =  1  in  f  expression,  the 

error  dynamics  (5.11)  and  the  control  (5.13)  transformed  into  the  error  dynamics  and  control 
obtained  from  regular  dynamic  inversion  formulation  (Padhi  et  al.  [2007]).  The  regular  dyna¬ 
mic  inversion  control  technique,  therefore,  is  a  subset  of  proposed  technique.  Furthermore, 
B^^KijI f  matrix  is  always  a  positive  definite  matrix,  hence,  the  nonlinear  error  dynamics  (5.11) 
is  always  be  stable.  Hence,  the  proposed  technique  is  named  as  state-constrained  generalized 
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nonlinear  dynamic  inversion  (SC-GDI)  control  technique. 

5.1.2  Neuro- adaptive  SC-GDI  Control  Formulation 

In  this  section,  The  formulation  of  adaptive  control  using  neuro -adaptive  SC-GDI  technique 
is  presented  for  an  uncertain  system.  The  desired  nominal  plant  and  approximated  model 
plant  dynamics  is  shown  in  equations  4.1  and  4.3  respectively.  Hence,  by  using  equation  5.11, 
the  nonlinear  error  dynamics  can  be  written  as 

Xa-Xa  +  Bl^KbifiXa  -Xd)  =  0  (5.14) 


here,  is  the  positive  definite  matrix.  Further,  equation  5.14  is  rearranged  in  a 

general  form  as 


f[X,Xa.Xd.U,UdJ[X))  =  0 


(5.15) 


here,  the  equation  5.15  is  the  general  expression  to  obtain  the  neuro -adaptive  SC-GDI 
control.  Further,  Ud  is  the  desired  control  formulated  in  subsection  5.1.1.  Moreover,  In  the 
expression  (5.15),  the  adaptive  control  U  is  obtained  by  using  Newton-Raphson  method 
provided  d{X)  is  known.  Further,  d{X)  is  obtained  by  weight  update  rule  explained  in  section 

4.2. 


5.2  Simulation  Results 

In  this  section,  the  simulation  results  are  presented  with  the  application  of  neuro -adaptive 
SC-GDI  control  on  the  actual  system.  The  model  uncertainty  is  assumed  in  the  aerodynamic 
forces  and  moments  coefficients  and  in  the  inertia  values.  The  disturbance  is  added  in  P, 
Q,  and  R  state  equations.  The  initial  condition  values  assumed  in  the  simulations  are  given 
as,  Vt  =  1789  m/Sy  a  =  2.94  degy  ^  =  0  degy  P  =  0  degIseCy  Q  =  0  degIseCy  R  =  0  degIseCy 
cf)  =  0  degy  6  =  2.94  degy  if/  =  0  degy  x  =  0m,  y  =  0m,  =  24000  m  and  initial  mass  Mi;  = 

136077.6  kg.  The  inner  loop  gains  fcp,  /cq,  and  kn  are  selected  as  2.0.  Moreover,  A  and  ks  values 
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(5.11)  for  P,  Q,  R  states  are  selected  as,  [0.001°,  0.001°,  0.001°]  and  0.1,  0.3,  0.1  respectively. 
Further,  the  intermediate  loop  gains  /c^,  kg,  and  are  selected  as  0.1.  The  outer  loop  control 
gain  kh  is  selected  as  0.01  and  further,  the  velocity  loop  controller  gain  is  chosen  as,  0.1. 


The  positive  definite  matrix  k^d  used  in  equation  4.4  is  selected  as  diag{5y  7,  5)  and 


kap  used.  The  disturbance  learning  rate  jt  is  selected  as  40  and  the  constraint  on  the  error 
(equation  4.6)  is  selected  as  O.ldegl sec.  Further,  the  simulation  is  carried  out  far  a  random 
disturbance  bounded  within  30%  of  model  uncertainty.  Moreover,  the  basis  functions  cf){x)  is 
selected  as  [1,  P,  P^,  Q,  Q^,  (PQ)^,  P,  P^,  (PP)^,  (QP)^,  (PQP)^]  and  the  lacobian  matrix  of 


the  basis  function 


is  already  given  in  equation  4.14. 


In  the  simulations,  three  type  of  plots  are  presented,  the  first  plot  is  obtained  by  the 
application  of  nominal  controller  on  nominal  dynamics,  further,  the  second  plot  is  obtained 
by  the  application  of  nominal  controller  on  actual  system  dynamics  and  finally,  third  plot  is 
obtained  by  the  application  of  adaptive  controller  on  the  actual  system.  It  can  be  seen  from  the 
plots  that  nominal  controller  is  not  effective  for  the  actual  system  and  results  into  an  unstable 
AHV.  The  velocity,  angle  of  attack  and  sideslip  angle  profiles  are  shown  in  Fig.  5.1  respectively 
with  the  application  of  adaptive  controller.  The  effectiveness  of  the  adaptive  controller  is  seen 
by  the  plots  where  the  desired  states  of  system  are  perfectly  tracked  by  the  actual  states  of 
system.  Further,  the  angle  of  attack  profile  is  bounded  within  the  limits  [-5°,  5°]  and  constant 
velocity  is  maintained  by  the  velocity  controller  to  prevent  the  inlet  unstart. 


Further,  the  altitude  of  the  AHV  is  maintained  close  to  the  reference  altitude  by  the 
application  of  SC-GDl  adaptive  controller,  as  shown  in  the  height  plot  of  Fig.  5.2.  The  mass 
profile  is  smooth  and  continuous  which  means  that  there  is  no  sudden  trust  requirement 
during  the  flight.  Moreover,  the  flight  path  angle  profile  is  stabilized  at  zero  by  the  application 
of  adaptive  control.  Further,  the  roll,  pitch,  and  yaw  angle  profiles  are  presented  in  Fig.  5.3 
respectively.  All  the  angle  profiles  are  bounded  and  the  error  between  the  desired  states  and 
the  actual  states  become  zero  in  finite  time. 
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Figure  5.1  -  Mach  Number  (M),  Angle  of  Attack  (a)  and  Sideslip  Angle  (yS)  Profiles  of  AHV 
using  Neuro-adaptive  Generalized  Dynamic  Inversion 


Further,  the  roll,  pitch  and  yaw  rate  profiles  are  shown  in  Fig.  5.4  respectively.  All 
the  angular  rate  profiles  are  converged  quickly  towards  the  desired  states,  which  means 
that  the  AHV  is  stabilized  by  using  the  derived  adaptive  controller.  Hence,  the  formulated 
neuro-adaptive  SC-GDI  controller  is  effective  for  controlling  the  uncertain  AHV  Further,  the 
adaptive  control  profiles  are  shown  in  Fig.  5.5.  The  control  surface  deflections  are  smooth  and 
continuous  (shown  in  zoomed  in  box)  and  within  the  saturation  limits  hence,  the  solution 
provided  by  the  SC-GDI  controller  is  physically  achievable. 
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Figure  5.2  -  Altitude  (/z),  Mass  (Mp)  and  Flight  Path  Angle  (y)  Profiles  of  AHV  using  Neuro- 
adaptive  Generalized  Dynamic  Inversion 
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Figure  5.3  -  Roll  Angle  (0),  Pitch  Angle  (0),  and  Yaw  Angle  (y/)  Profiles  of  AHV  using  Neuro- 
adaptive  Generalized  Dynamic  Inversion 
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Figure  5.4 -Roll  Rate  (P),  Pitch  Rate  (Q),  andYawRate  (R)  Profiles  of  AHV using  Neuro-adaptive 
Generalized  Dynamic  Inversion 
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Figure  5.5  -  Control  Profiles  obtained  using  Neuro-adaptive  Generalized  Dynamic  Inversion 
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Further,  the  simulation  results  are  presented  for  random  bounded  disturbances  between 
[-30%,  30%].  The  velocity,  angle  of  attack  and  sideslip  angle  profiles  are  shown  in  Fig.  5.6 
for  uncertain  AHV.  The  inlet  unstart  condition  is  successfully  prevented  by  the  proposed 
controller  because  the  angle  of  attack  is  perfectly  bounded  between  [-5°,5°]  and  velocity 
level  is  maintained  accordingly.  The  altitude  and  mass  profiles  for  AHV  are  shown  in  Fig. 
5.7  respectively.  The  smooth  and  continuous  profiles  are  obtained  in  presence  of  random 
disturbances  shows  the  applicability  of  the  adaptive  controller. 


Time  [sec] 


Figure  5.6  -  Mach  Number  (M),  Angle  of  Attack  (a)  and  Sideslip  Angle  (j6)  Profiles  for  Random 
Disturbances  using  Neuro- adaptive  Generalized  Dynamic  Inversion 


The  angular  angle  and  rate  profiles  are  shown  in  Fig.  5.8  and  Fig.  5.9  respectively. 
Minimal  transients  can  be  seen  in  all  the  plots  because  of  Jacobian  matrix  learning.  All  the 
angular  angle  and  rate  profiles  are  conversed  into  desired  trajectories  hence,  the  adaptive 
controller  performance  is  maintained  in  presence  of  AHV  model  uncertainties.  Therefore,  the 
high  fidelity  model  of  AHV  is  not  required.  Finally,  the  adaptive  control  profiles  are  shown  in 
Fig.  5.10,  where  all  the  trajectories  are  smooth  and  achievable. 
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Figure  5.7  -  Altitude  (/z),  Mass  (M^;)  and  Flight  Path  Angle  (y)  Profiles  for  Random  Disturbances 
using  Neuro-adaptive  Generalized  Dynamic  Inversion 
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Figure  5.8  -  Roll  Angle  (0),  Pitch  Angle  (0),  and  Yaw  Angle  iy/)  Profiles  for  Random  Disturbances 
using  Neuro-adaptive  Generalized  Dynamic  Inversion 
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Figure  5.9  -  Roll  Rate  (P),  Pitch  Rate  (Q),  and  Yaw  Rate  (P)  Profiles  for  Random  Disturbances 
using  Neuro-adaptive  Generalized  Dynamic  Inversion 
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Figure  5.10  -  Control  Profiles  obtained  for  Random  Disturbances  using  Neuro-adaptive  Gene¬ 
ralized  Dynamic  Inversion 
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5.3  Quantitative  Comparison  Between  Neuro-adaptive  DI  and  GDI 

In  this  section,  the  Monte-Carlo  simulation  results  obtained  from  the  application  of  neuro- 
adaptive  DI  and  neuro-adaptive  SC-GDI  on  the  actual  system  is  presented.  To  obtain  the  fair 
comparison,  results  are  simulated  for  the  same  disturbances.  The  Monte-Carlo  simulation 
results  are  shown  in  Table.  5.1,  where,  up  to  30%  disturbance,  both  control  techniques  works 
perfectly.  However,  as  the  disturbance  increases,  the  SC-GDI  proved  slightly  better  than 
the  DI  because  in  SC-GDI  technique  the  state  constraint  philosophy  is  used  along  with  the 
constrained  learning. 

Table  5.1  -  Quantitative  Comparison  Between  Neuro-adaptive  DI  and  GDI 


Model  Un¬ 
certainties 

Total  Trials 

SC-GDI: 

Successful 

Trials 

Percentage 
of  Success 

in  SC-GDI 

DI:  Success¬ 
ful  Trials 

Percentage 
of  Success 

in  DI 

10% 

3697 

3697 

100% 

3697 

100% 

20% 

3660 

3660 

100% 

3660 

100% 

30% 

3602 

3594 

99.78  % 

3580 

99.39  % 

40% 

3873 

3455 

89.2  % 

3371 

87.03  % 
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CONCLUSION 


Two  robust  neuro- adaptive  control  design  techniques  are  proposed  in  this  report  for  the  AHV 
in  the  cruise  phase  of  flight.  Further,  an  uncertain  nonlinear  winged-cone  model  is  used  to 
show  the  effectiveness  of  the  derived  controllers.  First,  DI  controller  is  used  to  obtain  the 
nominal  control  for  the  nominal  system  (no  disturbance).  The  effectiveness  of  the  nominal 
controller  is  shown  by  perturbed  initial  conditions. 

Further,  in  presence  of  model  uncertainties  and  disturbances,  the  robust  neuro -adaptive 
DI  controller  is  obtained  by  augmenting  DI  controller  with  proposed  constraint  disturbance 
learning  algorithm.  The  applicability  of  the  proposed  control  technique  is  shown  by  simulation 
results.  Further,  Monte-Carlo  simulation  results  are  also  presented  in  the  report. 

A  newly  developed  neuro-adaptive  SC-GDI  control  technique  is  also  proposed  in  this 
report  for  the  robust  control  of  uncertain  AHV.  BLF  is  used  in  the  formulation  of  SC-GDI  to 
obtain  the  state-constrained  control.  Further,  SC-GDI  control  is  augmented  with  the  proposed 
constraint  disturbance  learning  algorithm  to  formulate  the  neuro-adaptive  SC-GDI  control. 
Further,  the  effectiveness  of  the  proposed  neuro-adaptive  SC-GDI  control  is  validated  by 
Monte-Carlo  simulations. 

Finally,  It  has  been  shown  that  both  control  techniques  worked  perfectly  in  presence 
of  up  to  30%  model  uncertainties  and  prevented  inlet  unstart  along  with  the  stabilization  of 
the  AHV  Hence,  both  neuro-adaptive  DI  and  SC-GDI  can  be  used  to  control  the  air-breathing 
hypersonic  vehicle.  Moreover,  the  high  fidelity  model  of  AHV  is  also  not  required. 
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